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We present a mechanism for emission of electromagnetic terahertz waves by simulation. High 
T c superconductors form naturally stacked Josephson junctions. When an external current and a 
magnetic field are applied to the sample, fluxon flow induces voltage. The voltage creates oscillating 
current through the Josephson effect and the current excites the Josephson plasma. The sample 
works as a cavity, and the input energy is stored in a form of standing wave of the Josephson plasma. 
A part of the energy is emitted as terahertz waves. 

PACS numbers: 74.50.+r, 74.25.Gz, 85.25.Cp 



Continuous coherent terahertz waves have various ap- 
plications in scientific field such as biology and informa- 
tion science. One of the hurdles for technological ad- 
vancements in the terahertz region of electromagnetic 
wave is the development of sources for intense and con- 
tinuous coherent terahertz waves. Therefore, we inves- 
tigate a new mechanism for emitting intense continuous 
and frequency tunable terahertz waves. In the high tem- 
perature superconductors, the strongly superconduct- 
ing Cu02 layers and insulating layers are alternatively 
stacked along the c-axis of the crystals and form a natu- 
rally multi-connected Josephson junction called intrinsic 
Josephson junction (I J J). In the IJJ there appears a new 
excitation wave called Josephson plasma, the frequency 
of which is in the range of terahertai*2,. The frequency 
appears in the region inside the superconducting energy 
gap and the Landau damping is very weak, and thus the 
excited plasma decays by emitting a terahertz electro- 
magnetic wave. 

For investigating an emission mechanism of terahertz 
electromagnetic wave from the IJJ, we use the following 
model shown by Figure 1. In Fig. 1 the IJJ is shown in 
green and the electrodes of a normal metal (for example 
gold) are shown in yellow. An external magnetic field 
B applied in the direction of the y-axis induces fluxons 
in the direction. The centers of fluxons are in the insu- 
lating layers. In this system, the superconducting and 
normal currents almost uniformly flow in the direction 
indicated by J in Fig.l. The fluxons flow in the direction 
of the x-axis with a velocity v and induce the flow volt- 
age in the direction of the z-axis. These voltages creates 
the oscillating Josephson current along the z-axis by the 
Josephson effect, when temperature is low enough below 
T c and the superconducting current is smaller than the 
superconducting depairing current along the c-axis. This 
oscillating current interacts strongly with the Josephson 
plasma due to the nonlinear nature of the system and 
intensively excites the Josephson plasma wave as shown 
later. We use Bi 2 Sr 2 CaCu208+5 that is appropriate in 
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FIG. 1: Schematic diagram of a prototype model for ter- 
ahertz emission. Bi2Sr2CaCu20g+5 forms the IJJ shown in 
green, and electrode gold plates form the top and bottom 
electrodes shown in yellow. A dielectric waveguide shown in 
orange extends from the left surface of the IJJ. 



the experiments, and apply a magnetic field and external 
currents around J c the critical current to the IJJ. Then, 
the frequency of the plasma waves appears in the tera- 
hertz frequency range. The plasma wave is converted to 
an intense terahertz electromagnetic wave in the waveg- 
uide (dielectric) shown in orange in Fig. 1. 

In accordance with the mechanism mentioned above, 
we now derive the equations for the simulation. The 
superconducting order parameter of the Zth layer is ex- 
pressed as ipi(r,t) — A;e lw ( r '*). Since superconductivity 
in the Cu02 layers is strong at low temperatures, we as- 
sume A/(r, t) to be spatially and time independent. In 
this system, the current along the z-axis perpendicular 
to the Cu02 plane is given by a sum of the Josephson, 
quasi-particle, and displacement currents by 



Jz,i+i,i (r, t) = J c sin ipi+ij (r, t) + a c E 
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where J c is the critical current in zero magnetic field, a c 
is the normal conductivity along the c-axis, and ^ z .i+i,i is 
the electric field between (Z+l)th and Ith layers along the 
z-axis. The ipi+i t i is the gauge invariant phase difference 
defined as 
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with the vector potential A z (r, z, t) and the flux unit 0o- 
For the superconducting current densities in the Cu02 
plane, we use the generalized London equations, since the 
Ginzburg-Landau parameter is very large in cuprate su- 
perconductors. We insert Eqs. (f ) and (2) into Maxwell's 
equations taking account of the superconducting current 
in the afo-planes and the charging effect in the CuC>2 
plane. Following the calculation procedure^ and after 
cumbersome calculations, we have 

(1 - (A^)[{d^ l+ u + + sin *k+i,l) 

+a S i(d t ,(p' l+1 -p' l )+P(p' l+1 -p' l ))] 

= + dy,t/ji +lt i, (3) 

s'(l - aA^)p[ = dt'^i+^i - tPu-i). (4) 

In the above equation, we use normalized units for length, 
time, and charge density in the Ith CuC>2 plane, respec- 
tively defined as 

i x j i j / AgCJgp; 
x =— ,t =oj p t,andp l = : , (5) 

Jc 

where p\ is the charge density in the Ith Cu02 layer. 
The plasma gap frequency is given by lu p = c/( v / £c'A c ), 
e c being the dielectric constant along the z-axis and A c 
being the magnetic field penetration depth from the be 
surface plane. The parameters in Eqs. (3) and (4) are 
defined as 

, K b £cH 2 , s 4vrcr c A c 

C = TT ' a = ~Zf' s = t-, a^c? /? = -JTT (6) 
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where the s and d are the thickness of the superconduct- 
ing and insulating layers, respectively. X a b is the London 
penetration depth of the superconducting layer, and p is 
the Debye screening length. The operator A^ 2 ) is defined 
as A( 2 )/ ; = f l+1 - 2f t + fi—\. We use Eqs. (3) and (4) 
in the IJJ and use Maxwell's equations in the dielectric 
to simulate the emission of terahertz wave. 

If Eqs. (3) and (4) are linearized, we obtain the Joseph- 
son plasma solution. The Josephson plasma is a compos- 
ite wave of the Josephson current and electromagnetic 



wave. The amplitudes of the Josephson current and the 
electric field are always parallel to the c-axis. Depending 
on the wave propagation directions parallel to the c-axis 
and the a-axis, there appear respectively the longitudinal 
and transverse plasma waves. The transverse and longi- 
tudinal waves have been observed by Tamasaku et al^, 
and Matsuda et al(&, respectively. 

Let us show how the plasma wave is converted to the 
electromagnetic wave at the interface between the IJJ 
and the dielectric. As mentioned before, the plasma 
waves have the transverse and longitudinal components. 
However, the electromagnetic wave has only the trans- 
verse component. Therefore, only the transverse plasma 
wave can convert into the electromagnetic wave at the in- 
terface. When we solve Eqs. (3) and (4), we impose the 
following boundary condition. To connect the Josephson 
plasma wave in the IJJ to the electromagnetic wave in the 
dielectric at the interface, we put the usual electromag- 
netic boundary condition; the electric and magnetic fields 
parallel to the interface are continuous at the interface. 
The electromagnetic wave in the dielectric is assumed to 
transmit freely to outer space at the end surface of the 
dielectric. One possible way to achieve this condition is 
to make a gradation of the dielectric constant near the 
surface. It is assumed that the opposite surface of the 
IJJ is exposed to vacuum. 

Keeping in mind that the IJJ is f^S^CaC^Os+a and 
the dielectric is MgO, we chose A Q i,=0.4//m, A c =200yum, 
s=3l, d=12l, n=0.6A, a=0.l, /3=0.01~0.05, and the 
number of layers=20~1000, and take the dielectric con- 
stants along the z-axis in the IJJ and that of the dielec- 
tric wave-guide to be £=10. We apply a magnetic field of 
lTesla along the y-axis. We supply an external current 
from the gold electrodes shown in yellow in Fig. 1. We 
change the normalized external current J/J c from 0.2 to 
1.5 in step of 0.0125. The superconductivity in the top 
and bottom electrodes is assumed to penetrate to 0.075 
fim. We note that for the external currents mentioned 
above, the superconducting current part of J is always 
less than J c . The length of the IJJ is taken to be 100/im 
along the x-axis and the length of the dielectric is taken 
to be 50/mi along the x-axis. For each external current, 
the time evolution is simulated until the system reaches 
a stationary state, that is normalized time t'=600, and 
the emission power is calculated after that time. 

We assume that the system is uniform along the y- 
axis and make two-dimensional calculation in the x-z 
plane. We use the finite difference method to perform 
the numerical simulation. The simulation uses very large 
sized nonlinear equations heretofore difficult to compute; 
for a simulation using 10 6 spatial cells in the x-z two- 
dimensional model, it would take two years to simulate 
one case of 10 s time-steps using a personal computer with 
a 2GHz processor. Therefore, we used the Earth Simula- 
tor that has a peak performance of 40 teraflops and we 
carried out this simulation in one day. We simulated a 
total of sixty cases. Therefore, in sixty days we carried 
out simulations that would have required one hundred 
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FIG. 2: A snapshot of moving fluxons and the oscillating 
electric fields of the standing plasma wave. 
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and twenty years using a personal computer. 

Since the electromagnetic wave is the transverse wave, 
only the node-less Josephson plasma wave along the z- 
axis (the transverse plasma) can be converted into the 
electromagnetic wave with the same frequency. The max- 
imum stacking number that fulfill the condition is deter- 
mined by the applied magnetic field, and it increases as 
the magnetic field and A a & increase. We found that the 
stacking number that fulfill the condition is 100 layers 
under the magnetic field of 1 Tesla and X a b of 0.4/zm. 

Figure 2 shows a snapshot of the moving fluxons and 
the oscillating electric field. In the figure the static and 
uniform electric field have been subtracted. The black 
circles show the centers of the fluxons. The fluxons are 
almost randomly distributed and move uniformly in the 
x-axis direction. The speed of the fluxons is about 4% of 
the speed of light in the current J/ J c =0.6, /3=0.02 and 
under the applied magnetic field of 1 Tesla. The magnetic 
field of a fluxon spreads in the range of A c along the ab- 
plane. In the magnetic field of 1 Tesla, 120 fluxons are 
contained in the range. Therefore, the magnetic field of 
the moving fluxons are uniform and thus the fluxon flow 
create a spatially almost uniform voltage in the direc- 
tion of the z-axis and the voltage induces the oscillating 
current due to the Josephson effect and the oscillating 
current excites the Josephson plasma. The plasma wave 
is reflected at the interface between the IJJ and the di- 
electric, and at the surface exposed to vacuum, and it 
forms a standing wave in the stationary state as shown 
in Fig. 2. In the figure red and yellow show the electric 
field directing to the z-axis direction, and blue shows the 
electric field directing to the -z axis direction. 

Figure 3 shows a snapshot of the electric and magnetic 
fields of the standing wave at the 50th insulating layer 
from the bottom surface in the direction z-axis. The 
standing wave of the electric field along the z-axis oscil- 
lates with 2.86 THz around a constant electric field for 
the normalized current current J/ J c =0.8 and for /3=0.02. 



FIG. 3: A snapshot of the plasma wave that is a standing wave 
in the IJJ and the electromagnetic wave that is a propagating 
wave in the dielectric. 



The red line shows the interface between the dielectric 
and the IJJ. The components of the electric and mag- 
netic fields parallel to the interface fulfill the electromag- 
netic continuity condition. In Fig. 3 the green curves 
indicate the magnetic and electric fields at the time of 
1.2769ns after the system attains to the stationary state 
and the brown curves indicate the magnetic and electric 
fields at the time of 1.2772ns. The belly of an electric 
field oscillation with the largest amplitude is pinned at 
the interface, inducing the large amplitude of the electric 
and magnetic field oscillations in the dielectric. The ratio 
of the electric field to the magnetic field is equal to that 
of the electromagnetic wave at 4/im from the interface. 
The emission mechanism of the coherent terahertz elec- 
tromagnetic wave is similar to a laser mechanism. The 
IJJ itself works as a cavity. The energy is stored in the 
form of a standing Josephson plasma in the IJJ. A few 
percent of the energy is emitted as a terahertz electro- 
magnetic wave from the IJJ to the dielectric through the 
interface. 

In Fig. 4 we show the frequency spectra calculated by 
FFT analysis of the electromagnetic wave at a location 
4/im from the interface between the IJJ and the dielectric 
for the case of (3=0.03. When we change the normalized 
current J/J c from 0.85 to 1.3, the frequencies of the elec- 
tromagnetic waves change from 2.08 THz to 3.21 THz. 
The frequency spectra are very sharp. The frequency of 
the electromagnetic wave in the IJJ varies continuously 
by changing current J / J c , indicating that emission of the 
sample is frequency tunable. Some peaks have satellites 
due to the frequency modulation of the excited plasma 
wave. The frequencies calculated by the simulation are 
equal to the AC Josephson frequency estimated by the 
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FIG. 4: The frequency spectra analyzed by FFT for 10 nor- 
malized currents. The intensities are those of electromagnetic 
waves in the dielectric. The intensity is obtained by calculat- 
ing the Poynting vector at a location of 4/xm from the inter- 
face. 
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along the y-axis is 500/im. The emission occurs in some 
limited current region as seen in Fig5. When the current 
is smaller than a critical current, the emission power van- 
ishes. The reason is that nodes in the Josephson plasma 
appear along the z-axis, indicating inclusion of the lon- 
gitudinal component, and thus the emission power van- 
ishes as mentioned before. When the current is larger 
than the upper critical value, the fluxons take various 
dynamical orders 6 . In this state, the emission vanishes 
since the rigid fluxon lattice flow prohibits occurrence 
of the standing wave. The sharp dip of power appear- 
ing at J/J c =0.7 for (3=0.02 comes from the following 
reason. The fluxons flow with a small amplitude vibra- 
tion. When the Josephson frequency resonates with the 
fluxon vibration frequency and this resonance stabilizes 
the Josephson plasma, it depresses the conversion rate 
from the Josephson plasma to the electromagnetic wave. 
When (3 decreases, for example (3=0.01, the oscillations 
of electric and magnetic fields lose the coherence and the 
emission power vanishes. When (3 increases, the dumping 
becomes larger and the emission power decreases. 

In the present paper we performed the simulation of 
terahertz wave emission with the parameters correspond- 
ing to Bi2Sr2CaCu208+<5, and obtained the following re- 
sults: The sample works as a cavity, and the input energy 
is stored in a form of standing wave of the Josephson 
plasma. A part of the energy is emitted as terahertz 
waves. The emission of electromagnetic waves is in the 
form of continuous coherent terahertz waves with mW 
power in the present model. The frequency is tunable by 
changing the applied current. 



FIG. 5: The emission power (mW) of the electromagnetic 
wave as a function of the external current. 
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fluxon flow voltage within the accuracy of 1.3%. This 
fact suggests that the plasma wave is excided by the os- 
cillating current through the AC Josephson effect. 

Figure 5 shows the emission power (the Poynting vec- 
tor) measured at the location 4/im from the interface as a 
function of the external current. The length of the sample 
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